Frustrated Bose-Einstein condensates with non-collinear orbital ordering 



(N 

o 

(N 
< 



<Z3 

^— > ' 

c : 



c 

o 

> 

m 
O 

in 

O 



X 



Zi Cai/'B Yu Wang, 2 '0 and Congjun Wu 1 - 2 

'Department of Physics, University of California, San Diego, CA92093 
2 School of Physics and Technology, Wuhan University, Wuhan 430072, China 

We investigate the unconventional Bose-Einstein condensations with the orbital degree of freedom 
in the 3D cubic optical lattice, which gives rise to various exotic features absent in conventional scalar 
and spinor Bose-Einstein condensations. Orbital angular momentum moments are formed on lattice 
sites breaking time-reversal symmetry spontaneously. Furthermore, they exhibit orbital frustrations 
and develop a chiral ordering selected by the "order-from-disorder" mechanism. 
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Bose-Einstein condensation (BEC) is a striking phe- 
nomenon of quantum many-body systems, character- 
ized by a uniform phase that spontaneously breaks the 
U{\) symmetry. By introducing extra degrees of free- 
dom, novel quantum condensates with even more exotic 
symmetry breaking patterns and topological structures 
emerge. A familiar example is the superfluid 3 He phases, 
which is a spin-triplet p-wave Cooper pairing conden- 
sate with both spin and orbital degrees of freedom 1 * 2 -. It 
exhibits a variety of rich structures that simultaneously 
incorporate the symmetries of liquid crystals, magnets 
and scalar superfluids. Consequently, the superfluid 3 He 
systems possess fundamental connections with particle 
physics, and exemplifies fundamental concepts of mod- 
ern theoretical physics 3 . 

The rapid developments of cold atom gases provide an 
opportunity other than 3 He to explore exotic condensa- 
tions with internal degrees of freedom. Spinor atomic 
gases, composed of atoms with hyperfine spins, simul- 
taneously exhibit magnetism and superfluidity 4-6 . Fur- 
thermore, orbital is a degree of freedom independent of 
spin and charge. It is originally investigated in condensed 
matter transition metal oxides, which plays an important 
role in superconductivity, metal-insulator transition, and 
quantum magnetism 7,8 . Introducing orbital into cold 
atom gases has been theoretically investigated, which 
leads to unconventional BECs with complex-valued con- 
densed wavefunction of bosons and spontaneously break- 
ing of time-reversal symmetry*^—. Excitingly, the recent 
experimental progress has realized the meta-stable BECs 
in high-orbital bands exhibiting complex- valued conden- 
sate wavcfunctions at non-zero wavevectors and orbital 
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In this paper, we investigate the properties of an or- 
bital BEC in the p-orbital bands of a cubic optical lattice. 
Although orbital BECs share many properties with the 
ferromagnetic phase in spinor BECs, there are crucial 
differences between them. In spinor BECs, the internal 
degrees of freedom of hyperfine spin degree is not cou- 
pled to the lattice. However for orbital BECs, the or- 
dering of the orbital angular momentum comes from the 
atom motion within each optical site, and thus is closely 
related to the motion of atoms in optical lattice. As we 
will show, the uniqueness of the orbital degree of freedom 



gives rise to a whole host of exotic phenomena, such as 
orbital frustration and concomitant non-collinear orbital 
orderings selected by the "order-from-disorder" mecha- 
nism. The selected ordering pattern exhibits an orbital 
angular momentum moment chirality. 

The Hamiltonian of bosons pumped to the p-orbital 
bands of a cubic optical lattice is described by a multi- 
orbital Bose-Hubbard model, H = H t + H int , 
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where /i, v — x,y, z denote the orbital indices, a is the 
lattice constant, p^^ ( p* n S) are annihilation (creation) 
operators for bosons at site f in orbital fi. n? is the to- 
tal particle number operator and L? represents the total 
orbital angular momentum on site r. fn and t± describe 
the nearest-neighbor hopping matrix elements along the 
longitudinal and transverse directions, respectively. Us- 
ing the terminology of chemistry, they are denoted as a 
and 7r-bonding, respectively. Due to the odd parity of p- 
orbitals, t» and t± are positive. The strong anisotropy of 
the p-band Wannier wave-function implies that tj_ -C h\ . 
The onsite interaction term, Hint, reflects the Hund's 
type physics generalized to bosons, i.e., bosons prefer 
to occupy complex-valued orbitals of the Pe t + ipg 2 tyP e 
with ei _L e.%. This complex- valued orbital has larger 
spatial extension than the real ones of p u , and thus the 
repulsive interactions are reduced and simultaneously the 
onsite orbital angular momenta are maximizec&i 3 -. 

We consider the orbital superfluid phase. A remark- 
able feature of the band structure is that the energy 
minima of p M (/i = x, y, z)-orbitals are located at finite 
momenta Q n rather than at zero momentum, which are 
Qx = (5,0,0); Q v = (0,f,0); Q z = (0,0,f) for the 
three p-orbital subbands, respectively. In the 3D cu- 
bic lattice, we will show that the onsite orbital angu- 
lar momenta are no longer collinear but exhibit orbital 
frustrations. The single-particle states ipQ = e lC "^' r 
(fj, — x, y, z) are degenerate, thus any condensate wave- 
function of a linear superposition of these states, 
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fh, n and I = rh x n, which is just the SO (3) group space 
and can be expressed in terms of Euler angles (0, 9, 7), as 
illustrated in Fig. HJa). Note that the multiplication of 
an overall U(l) phase e~ llp is equivalent to the rotation 
of the triad around / by the angle ip. Therefore, the U(l) 
superfluid phase is absorbed into the SO (3) group con- 
figuration space. For a given triad configuration mo, no 
and lo, the corresponding real space distribution of the 
OAM orientation L(r) becomes 



FIG. 1: (a) The S0(3) manifold of the 3D p-orbital BEC order 
parameter in terms of Euler angles, where 1 is the direction of 
the orbital angular momentum; (b) Sketch of the non-collinear 
orderings of orbital angular momenta in real space from Eq. 
©• 



yields the same kinetic energy. The complex vector c = 
(01,02,03) satisfies the normalization condition, |c| 2 = 1. 
Next we will consider the interaction effect to lift the 
degeneracy and select the condensate wavef unctions. 
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FIG. 2: (a) 2D orbital BEC in the square lattice without 
frustration, (b) A typical configuration of 3D orbital BEC in a 
cubic lattice with frustration. The thick blue bonds minimize 
both the transverse and parallel hopping energy, while the 
thin red bonds only minimize the transverse hopping energy. 

The SO (3) degeneracy at the classical level: At the 
classical level (neglecting quantum fluctuations so that 
the boson operator can be replaced by its average value) , 
Hi n t is minimized if the coefficient vector c in Eq. @ can 
be expressed as c = (fh + in) , where rh and h arc two 
mutual perpendicular unit vectors. Transforming back 
into real space, for the lattice site with integer- valued co- 
ordinates r = (r x ,r y ,r z ), its onsite orbital configuration 
is ^(p fe 



ipi 2 ) with the relation 



ex = {P x ) r *{Py) r *{P z ) r >m; 

£ 2 = (p x ) r °>(p y y»(p z ) r *fi, 
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where P x , y ,z are reflection operator with respect to x, y, z- 
axes, respectively, e^a remain orthogonal to each other, 
and the onsite orbital angular momentum L(r) \\ e\ x §2, 
such that L 2 is maximized to minimize Hi nt . This de- 
notes that at the classic level, the ground state manifold 
is just the configuration space of the 3D orthogonal triad 



which is non-collinear as shown in Fig. [lib). 

Similarly to the case of frustrated magnets, this classic 
level degeneracy is a consequence of orbital frustration, 
which means that it is impossible to find an orbital con- 
figuration that simultaneously minimizes the energy of all 
the bonds in the lattice. To illustrate this point, let us 
recall the previously studied 2D case for a comparison-. 
The staggered OAM configuration in Fig. [5] (a) simulta- 
neously minimizes both the parallel (in) and transverse 
(fx) hopping energies at all bonds of the square lattice, 
and thus there is no frustration. However, in the 3D cu- 
bic lattice, this is no longer the case. For example, if we 
take a similar state \Q xy ) = -^(\Qx) + AQy))i as shown 
in Fig. [UJb), the hopping energy of all bonds along x 
and y directions can be minimized, but the cr-bond along 
the z-direction is broken. Since the hopping Hamilto- 
nian Eq. [T]does not preserve the S*0(3) symmetry, thus 
this classic level degeneracy should be lifted by quantum 
fluctuations. 

Order- from- disorder In frustrated magnetism, the infi- 
nite degeneracy is usually lifted by quantum or thermal 
fluctuations, which is known as "order-from-disorder" 
mechanism22i2i. Below we perform the same analysis 
to the 3D p-orbital BECs. If we take quantum fluctu- 
ations around the mean-field values into account: — 
(p a ) +5p^ and calculate the fluctuation-corrected ground 
state energy, quantum fluctuations lift the S0(3) classical 
degeneracy. We consider two typical condensate config- 
urations and compare their ground state energies (the 
reason to choose these two states is due to their high 
symmetry) , 



\Q dlag ) = —(\Q x )+ e ^\Q y )+e-^\Q z )), (5) 
\Q Z ) = -L(\Q x )+i\Q y )). (6) 

In the state of \Qdiag), OAMs are along the body- 
diagonal directions, while for the state of \Q Z ), OAMs 
are along the z-direction. These two configurations are 
degenerate at the classic level. 

Here we perform the standard Bogoliubov analysis^ to 
calculate the zero-point motion energy of quasi-particles 
for these two configurations. We use the state of \Qdiag) 
as an example, and the calculation for \Q Z ) is rather sim- 
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FIG. 3: The energy difference between the orbital BEC with 
OAM towards z-direction and body-diagonal direction AE = 
E z — Ediag with t\\ =t,tj_ = 0.05t and the filling factor n=2. 



ilar. For each p-component the order parameter is 

<p B ) = (-1)^0, ipy) = {-lyve^cf 

(p z ) = (-iPe-^</, (7) 

To calculate the Bogoliubov spectra, we consider quan- 
tum fluctuations around the mean-field values: p M = 
(Pfj.) + fiPv Expanding to the quadratic level, we arrive 
at 



ih?^- = M(k)*(k), 



(8) 



where ^(k) is a 6-component vector that represents 
the fluctuations: *(k) = [Sip(k), S^(-k)] T , SifiQn) = 
[5p x (k),5p y (k + Q xy ),5p z (k + Q xz )}, Q xy = (n,ir, 0), 
and Q xz = (n, 0,n). M(k) is a 6 x 6 matrix as 



U(k) A(k) 
-At(k) -K(-k) 
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in which both H and A(k) are 3x3 matrices: 
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where u; = git^; e£ = 2 I] l /[ t || < W~ < -L( 1 ~<W] cos(fc„a) 
is the single particle energy spectrum for the \i band bo- 
son. The self-consistent equation to determine value of 
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where e(k) = (eg + 4+Q xy + e£ +Qx J/3, and n is the fil- 
ing factor. The contribution from the zero point motion 
energy to the ground state energy can be written as 



E° diag = -ZUnl -Un c -t+lj2 \A"(k) (e(k) + 4[/n c 



where n c — |0d 9 | 2 , t = t\\ + 2t±. Performing the same 
process, we obtain the correction for \Q Z ), and the dif- 
ference AE — E z — Ediag is plotted in Fig. [3] 

For fixed parameters U, tu , tj_ and boson density n, the 
energy of \Qdiag) is always lower than that of \Q Z ) i.e., 
AE = E z — Ediag > 0, which means that orbital BECs 
in a cubic lattice prefer to develop OAM moments along 
the body-diagonal directions. Such a configuration has a 
high symmetry, and all the bonding strengths are uniform 
in the lattice. In comparison, all the er-bonds along the 
z-direction are broken in the state of \Q Z ). This "order- 
from-disorder" phenomenon is another feature that dis- 
tinguishes orbital BECs from spinor BECs. 

Spontaneous chiral orbital order: In condensed mat- 
ter physics, the presence of spin chirality plays im- 
portant roles in frustrated magnetism- 3 , doped Mott- 
insulators^, and the anomalous quantum Hall effect 25 . 
Here, we find that the orbital BECs in Eq. ([5]) spon- 
taneously develops a chiral orbital angular momentum 
ordering. Eq. ([5]) has a time-reversal partner as 



\Q' dl ag) = ^=(\Q, 



)+e- l ^\Q y )+e^\Q z )). (11) 



The orbital angular momentum orderings of Eq. [5] and 
Eq. [IT] are plotted in Fig. HI T° distinguish the chirality 
between \ Qdiag) and \Q' d i ag ), we define a nonzero chirality 
(similar to the case of chiral spin liquid 2 ^) : 



Xijk = k ■ \lj X Ik), 



(12) 



where ijk denotes three sites of the four corners of the 
plaquette in a clockwise direction, as shown in Fig.UJ(c). 
In contrast to the chiral spin liquid 2 ^ (where (S) — 0, \ 7^ 
0) and the spinor BEC (where (S) 7^ 0, \ — 0), the orbital 
BEC in our case simultaneously exhibits non-vanishing 
orbital order and orbital chirality order ((I) ^0,x^0). 
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FIG. 4: Ordering patterns in a plaquette for (a) \Q\ iag ) and 
(b) \Qdiag) with opposite chirality; (c) Definition of the chi- 
rality in a plaquette. 

Excitations: Now we discuss the collective modes and 
elementary excitations of the orbital BECs in the state 
of \Qdiag)- From the EOM Eq. (JSJ, we obtain the 
collective modes by diagonalizing the matrix A4. In 
the long- wavelength limit (k — > 0), we find that there 
are three modes: (1) the gapless Goldstone mode cor- 
responding to the fluctuation of the superfluid phase 
of the orbital BECs with linear dispersion, o>i(k) sa 
y/2U\(t)\ 2 {2tj_ + t||)|k|, (2) the orbital wave mode cor- 
responding to the fluctuation of the OAM around its 
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ground state directions, W2(k) w — -k 2 , (3) the 
gapped mode, which describes orbital excitations cor- 
responding to the flipping of orbital angular momenta, 
w 3 (k) « ^tifc 2 + 2[/|0| 2 . 

Ginzburg- Landau (GL) free energy: To better under- 
stand the various phases in orbital BECs, we first identify 
the order parameters of the different phases. The orbital 
BECs are characterized by both superfluidity with order 
parameters <p M (/i = x, y, z) and orbital order L, which is 
a 3-component real vector (after the transformation 
defined in Eq. |4j the non-collinear order can be trans- 
formed into a spatially uniform one). Notice that these 
two different orders are not necessarily simultaneously 
present. To clarify this point, we construct a GL free 
energy in terms of x M and L: 

? = £^l 2 + ^l 4 + r'W 2 + "'(|£| 2 ) 2 

the last term is the minimal coupling between the super- 
fluid and orbital orders and the parameters <?, v! , u are 
positive in our case. For the the orbital BEC in Eq. [5J L 
is along the body-diagonal direction and the free energy 
Eqll3l can be simplified to 

T = r a <p 2 + u </> 4 + r'J 2 + uX - 9o<p 2 l + • ■ • • (14) 

The free energy J- describes both thermodynamic phase 
transitions which are driven by temperature and quan- 
tum phase transitions where the parameters are a func- 
tion of U/t. By minimizing the free energy in Eq. (Q~4|), 

we find that: for r' Q < and ro < \J~^^i both (j> and 
I are non-zero at the free energy minima, which corre- 
sponds to orbital BECs with both superfluidity and or- 
bital order; for r' < and r > \J—^r-, <j> = while 
I 7^ 0, so in this case the single particle condensations are 
suppressed by thermal or quantum fluctuations while the 
orbital order is preserved; for r' > and ro > 0, both <p 
and I are zero, which means that both the superfluidity 
and orbital order have been destroyed, corresponding to 
the high temperature normal phases or featureless Mott 
insulator at zero temperature. 

Experimental detection: Next we discuss the experi- 
mental detection of orbital ordering and unconventional 



BEC characterized by the ordering parameters in Eq. . 
The condensate wavefunction in Eq. lfTTj) is a superposi- 
tion of the single-particle states of three condensate mo- 
menta with equal weights but different phases, as a con- 
sequence, the time-of-flight (TOF) image will exhibits 
three peaks with the same height in the points corre- 
sponding to Q x ,y,z- However, TOF images only pro- 
vide the single-particle density distribution in momentum 
space, while the key information about the relative phase 
between different condensate components is lost during 
TOF. To measure the phase difference e ~ defined in 
Eq.flT]), the phase sensitive detections proposed recently 
can be employed 2 ^. Impulsive Raman operations can be 
applied to couple any two among the three components 
at different momenta. The relative phase information can 
be read out from the interference pattern from the TOF 
imaging. 

At last, we will briefly discuss the orbital BECs in 
higher bands. Recently, an unconventional BEC in the 
/-band of a bipartite optical square lattice has been 
observed experimentally^. Surprisingly, d-band BECs 
have also been observed in a distinct field: the exciton- 
polariton condensate 2 ^. Orbital BECs in the p-band of 
3D optical lattice have three components (p X) p yi p z ), 
and the interactions favor a ferro-orbital state with OAM 
L = I, which makes it similar to the ferromagnetic phase 
of spinor BECs with F = 1. Analogously, orbital BECs 
in higher bands behave similarly to spinor BECs with 
higher spin2&~— . Apparently for orbital BECs in higher 
bands, the geometry and symmetry group of the order 
parameters is far more complex and may give rise to 
richer physics. 

In conclusion, we investigate the frustrated orbital or- 
dering of p-band unconventional BECs in the cubic lat- 
tice, and we find that the uniqueness of the orbital de- 
gree of freedom gives rise to a lot of interesting phenom- 
ena that absent in the spinor BECs and superfluid 3 He, 
such as orbital frustration and concomitant non-collinear 
orbital orderings selected by the "order-from-disorder" 
mechanism. The chiral symmetry breaking and the ele- 
mentary excitations have also been discussed. 
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